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We give a general review of recent developments in the theory of vortices in superfluids and 
superconductors, discussing why the dynamics of vortices is important, and why some key results 
are still controversial. We discuss work that we have done on the dynamics of quantized vortices in 
a superfluid. Despite the fact that this problem has been recognized as important for forty years, 
there is still a lot of controversy about the forces on and masses of quantized vortices. We think 
that one can get unambiguous answers by considering a broken symmetry state that consists of one 
O l_ vortex in an infinite ideal system. We argue for a Magnus force that is proportional to the superfluid 

density, and we find that the effective mass density of a vortex in a neutral superfluid is divergent 
at low frequencies. We have generalized some of the results for a neutral superfluid to a charged 
system. 

I. INTRODUCTION 

> 

^v^j ' From the. very beginning it has been realised that quantized vortices play an important part in the behavior of 
y—t [ superfluidsEJ. Both in neutral superfluids and in superconductors it is the vortices that provide a mechanism for 
ON 1 the decay of superfluid currents in a ring. The circulation, for a neutral superfluid, or the trapped flux, for a 
superconducting ring, is quantized, and the current can only decay by a change of the quantum number by an integer, 
which can occur by the passage of a vortex (or quantized flux line) across the ring, from one edge to the other. In 
superconductors in a high magnetic field, the motion of flux lines is the main mechanism for electrical resistance. At 
high temperatures the movement of vortices is a thermally activated process, but at low enough temperatures the 
dominant mechanism must be by quantum tunneling. It is therefore important to understand the dynamics of vortices, 
in order to be able to evaluate the dissipative processes that occur in neutral superfluids and in superconductors. 
■ Despite the obvious importance of the problem, the theory has been in a most unsatisfactory state. There are 
many conflicting results in the literature. We did not realise the extent of this disagreement, and it was initially a 
O ■ surprise to us that a result we obtained would be dismissed by one knowledgeable critic as too obvious to be worth 
discussing, and by another equally eminent critic as well known to be wrong; this has happened to us several times. 
^ , There are real problems here, connected with questions of suitable boundary conditions, and there is often a question 
of whether two different calculations are finding the same result by two different ways, or if they are finding two 
different contributions which must be added. 

To our surprise, we found, about two years ago, that we could get an exact result for one of the two parameters that 
determine the transverse force on a moving vortex, using only general properties of superfluid ordcrfl. The second 
parameter was determined by a straightforward thermodynamic argument a little over a year latera. The first of these 
results still seems to be controversial 4 ], and some elements of the argument deserve closer scrutiny than they have 
received so far, but it is our belief that it should be possible to construct a firmly founded theory on the basis that 
we have tried to establish. 
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II. ELECTRONS IN MAGNETIC FIELDS AND VORTICES 



There are strong analogies between the behavior of electrons in strong magnetic fields and of vortices in superfluids. 
These analogies enable us to make use of some of the insights that have been obtained in the study of the quantum 
Hall effect to understand problems connected with vortex dynamics. 

In both cases there is a transverse force proportional to velocity. The Lorentz force for electrons is proportional 
to the vector product of the electron velocity and the magnetic field, Fj, = — ev x B. This can be represented by 
a term eBxy in the Lagrangian. The Magnus force acting on a vortex is proportional to the vector product of the 
velocity of the vortex relative to the fluid and a vector directed along the vortex core. Each of these forces leads to 
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a path-dependent but speed independent term in the action. In a quantum theory the phase is equal to the action 
divided by h, this cospesponds to a Berry phase, a phase which depends on the path, but not on the rate at which 
the path is traversedtffl. 

In both cases there is considerable arbitrariness in the value of this phase. For the electron this is due to the 
arbitrariness of the vector potential which is used to represent the magnetic field, while for the vortex there is a 
similar arbitrariness in the way the transverse force is represented in a Lagrangian. In either case the change in action 
or phase when the electron or vortex traverses a closed path is well determined. For the electron the phase change on 
a closed path is equal to 2ir times the number of flux quanta enclosed by the path. For the vortex the phase change 
is equal to 2tt times the average number of atoms enclosed by the surface swept out by the closed path of the vortex. 

The dominant correction to quantization of the Hall conductance comes from tunneling or activated transport 
between states on the two edges of a quantum Hall bar. The dominant mechanism for decay of supercurrents is 
tunneling or activated transport of vortices across the system. In many real systems there is a tangle of pre-existing 
vortices frozen in when the system is cooled below the critical temperature, and these can serve as sources for the 
vortices that cross the system. Under ideal conditions, and some modern experiments on helium approach such ideal 
conditional, there are no vortices in equilibrium, and a vortex loop must be created from nothing in the interior, or 
a line must be created at the boundary (with its image this constitutes a loop), and cross the system to be annihilated 
at the opposite boundary. 

Electrons in a magnetic field have a fast cyclotron motion around the guiding center. Canonical variables can be 
chosen as the two pairs v x ,v y , rescaled by m 2 /eB, which give the fast motion, and the guiding center coordinates 
X,Y, rescaled by eB, which give the slow motion. Vortices also have such a fast cyclotron-like motion, in which the 
vortex core circles around the center of the flow it induces. In addition, since the vortex is a string rather than a 
point, it has the low frequency Thomson modes, circularly polarized modes of oscillation analogous to the modes of 
oscillation of a stretched string. For the vortex the two coordinates X, Y of the position of the vortex in the plane 
perpendicular to the core are also conjugate variables, as is manifest from the classical theory of vortex motion which 
can be found in Lamb's Hydrodynamics, and which Lamb credits to an 1880's book on Mechanics by Kirchhoff. 

The most important difference is that we think we understand the Schrodinger equation for electrons, whereas a 
vortex in a superfluid is a complicated many-body entity. The relation of the collective variables describing the vortex 
motion to the single-particle variables describing the superfluid is not obvious. 

III. GENERAL FEATURES OF VORTICES IN SUPERFLUIDS 

A vortex is a composite object in a many-body system. Its motion may be described by collective variables, but its 
structure depends on all the single-particle variables of the superfluid,_and the relation between these single-particle 
variables and the collective variables is, as usual, obscure. FeynmanEI proposed describing a vortex by taking the 
ground state wave function, symmetric in all the single particle variables for a boson superfluid, and multiplying it 
by a factor of the form 

JV 

where 8j is the azimuthal angle made by the particle j with the vortex core, rj is its distance from the core, and / 
is some real function which is close to unity everywhere except where rj is of the order of the radius of the vortex 
core, and which goes to zero at rj — in order to prevent large kinetic energy contributions from the rapid variation 
of phase close to the vortex core. A similar description, of the core is obtained for the Ginzburg-Pitaevskii equations 
for the order parameter near the critical temperatureEj, or fr om t he Gross-Pitaevskii nonlinear Schrodinger equation 
for the condensate of the dilute Bose gas at zero temperaturdlilO. In these theories the vortex is described by mean- 
field-like equations, so that the position of the singularity at the vortex core has a sharp value, although we know that 
the two components of its position in the plane are conjugate variables. Somehow we should be able to construct a 
quantized version of the theory taking account of this result of the Magnus force. 

In a strongly type II (Shubnikov) superconductor the situation is somewhat similar, except that the current circu- 
lating round the vortex core generates a magnetic field parallel to the core, which in turn generates a vector potential 
that reduces the current, so that a total of one quantum of flux h/2e is associated with the vortex, or flux line, and 
no current is associated with the change of the phase angle at large distances. In a type I (Pippard) superconductor 
the character of the singularity is mainly trapped flux, but the singly quantized flux line is not stable in a uniform 
magnetic field, and it is thermodynamically favorable for the flux lines to aggregate and form a region of normal 
metal. In either case Landau-Ginzburg theory can be used to describe the vortex. 
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In classical incompressible fluid mechanics the hydrodynamic mass of a vortex is of order of mass of fluid displaced, 
but it depends in detail on the core structure. Since vortices in low temperature superfluid helium are measured to 
have a rather small vortex core radius, smaller than the average interatomic spacing, this mass density is relatively 
small, and is taken to be zero in some-calculations. In recent work Duan and Leggett showed that the inertial mass of 
a vortex.in a superconductor is finite£3, but Duan argued that the mass density of a vortex in a neutral superfluid is 
infinitaij. He originally described this as a result of the quantum nature of the fluid, and we found this very hard to 
accept. Actually it is true of all compressible fluids, but the divergence is logarithmic in the frequency of the motion, 
with quite a small coefficient, as Demircan, Ao and Niu have pointed outtJ. Under realistic circumstances, such as 
in the free cyclotron motion of the vortex, or in vortex tunneling, the logarithm may be quite small, and this term 
relatively unimportant. 

In liquid helium at relatively high temperatures, close to the critical temperature, the largest force on a moving 
vortex, or on a vortex that is held still while the fluid streams past it, is likely to be a drag force due to the scattering 
by the vortex of the excitations that make up the normal fluid. At lower temperatures the transverse (Magnus) force 
should dominate, but the understanding of the Magnus force is complicated by the existence of the two components of 
the fluid, which may affect the vortex very differently. As we have discussed already, the Magnus force has important 
implications for the Berry phase, and for the quantum uncertainty of the position of the vortex. 

For a superconductor the situation is far more complicated, since not only is there the magnetic field due to the 
motion of the electrons to be considered, but the effects of disorder in the positive background are vital. Disorder makes 
the conductivity of the normal metal finite, and produces a drag force on vortices even at rather low temperatures, 
but also, if the disorder is on a large enough scale, pins the vortices and reduces the flux flow resistivity. 

In our work we have concentrated on understanding an isolated vortex in an ideal, uniform, infinite superfluid. 
Our aim has been to understand the parameters that come into the dynamics of a vortex when its velocity relative 
to the background fluid is small — the effective mass, the transverse component of the force, and the longitudinal 
(dissipative) component of the force. This is clearly not a program for a complete understanding of vortex dynamics, 
since, even if it were completely successful, we might still be concerned with strongly nonlinear regions in realistic 
situations, such as those found when quantum tunneling appears to be observed. Particularly for the transverse force, 
we think we have clean and precise results that are — inevitably — in conflict with widely accepted theories. 

In our work the vortex is controlled by some pinning potential that can be manipulated from outside. The pinning 
potential can be rather weak, or a macroscopic wire, so long as it has cylindrical symmetry. For quantities such as 
the effective mass and the longitudinal force the nature of the pinning potential has an effect on the answer, and we 
may need to consider some suitable limiting process to make the strength of the potential tend to zero, but for the 
transverse force we find that the answer is independent of the form or strength of the pinning potential. 

IV. THE MAGNUS FORCE IN NEUTRAL SUPERFLUIDS 

There is no agreement about what the forces acting on a vortex in a neutral superfluid are. The simplest quantity 
to calculate should be the component of force perpendicular to the motion of the vortex relative to the substrate, the 
analog of the Magnus force for classical fluids, yet two recently quoted forms look quite different. In Donnelly's book 
on Quantized Vortices in Helium I&Q he quotes the force per unit length as 

F f = K X [p s (Vy - V s ) + (fin - cr)(vy - V„)] , (2) 

where K is a vector along the vortex line whose magnitude is the quantum of circulation h/M, vy, v s and v n are the 
velocities of the vortex, the superfluid component and the normal fluid component, and p s and p n are the superfluid 
and normal fhiid-densities; a is a coefficient whose value is not exactly determined. Volovik, however, in a number of 
recent papersO'EJ, quotes the form 

F t = K x [p(v y - v 8 ) + p»(v s - v„) + C F (v„ - Vy)] , (3) 

where the term with coefficient Cf occurs only for fermion superfluids, such as the B phase of superfluid 3 He, and is 
due to spectral flow of the low energy states in the vortex core. The first term in each of these expressions is referred 
to as the Magnus force, the term proportional to p n as the Iordanskii force, so the use of these two terms is quite 
different for the two authors. Donnelly's term proportional to a comes from phonon or roton scattering by the vortex, 
and it is only if this is equal to zero that the two expressions are in agreement for the case of bosons. 

Whatever the form of this force, Galilean invariance tells us that there are only two parameters to be determined. 
If we know the coefficients of vy and v s , the coefficient of v n must be equal to minus the sum of the other two 
coefficients. We argue, in the rest of this section, that the only transverse force has the form 
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F t = p s K x (vy - v s ) , 



(4) 



by determining separately the coefficients of v s and vy. 

Wexlem has given a thermodynamic argument to show that the coefficient of K x v s is indeed — p s . This result 
seems to be uncontroversial, and is in agreement with both Donnelly and Volovik. The argument is essentially a 
thermodynamic argument, which considers a reversible change of the circulation in a ring by moving a vortex slowly 
across the system under equilibrium conditions. 



vortex 




FIG. 1. Vortex being moved from the outside to the inside of a ring of superfluid, increasing the superfluid circulation around 
the ring by one quantum unit. 

Consider a macroscopic ring, such as the one shown in fig. 1, with average radius R, width (difference between 
outer radius and inner radius) L y and height L z . For simplicity we assume L y << R, but this is not essential, and 
the result is independent of the shape of the ring. Initially there are n quanta of circulation trapped in the ring, 
giving superfluid velocity v s = nKo/27ri?, and the normal fluid velocity is zero, since the boundaries of the ring are 
stationary. A pinning potential is used to insert adiabatically one vortex, which is created at the outer boundary, 
moved slowly across the system under constant temperature conditions, and annihilated at the inner boundary. The 
effect of this extra vortex is to increase the circulation from n units to n + 1, increasing the superfluid velocity by 
Sv s — kq/2ttR. This increases the free energy by 

AF = (2irRL y L z )p s v s 5v s = p s K v s L y L z , (5) 

since superfluid density is defined in terms of the free energy change when the superfluid velocity is changed. This 
must be compared with the work done in moving the vortex of length L z a distance L y isothermally across the ring, 
which is 

W = F t L y L z . (6) 

Comparison of these two shows that the magnitude of the transverse force per unit length, under conditions in which 
v n and vy are both negligible, is 

\F t \ = PsKqv s . (7) 

More careful analysis gives the sign and direction of this force as 

Ft - - Ps K x v s . (8) 

This argument determines the coefficient of v s in the transverse force. 

To determine the coefficient of Vy, Thouless, Ao and Niutl consider an infinite system with superfluid and normal 
fluid asymptotically at rest (v„ = = v s ) in the presence of a single vortex which is constrained to move by moving 
the pinning potential. For simplicity we describe the two-dimensional problem of a vortex in a superfluid film, but the 
three-dimensional generalization is straightforward. Also we restrict this discussion to the ground state of the vortex, 
but the generalization to a thermal equilibrium state is straightforward. The reaction force on the pinning potential is 
calculated to lowest order in the vortex velocity vy . This can be studied as a time-dependent perturbation problem, 
but this can be transformed into a steady state problem, with the perturbation due to motion of the vortex written 
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as ivy ■ grad . The force in the y direction on a vortex moving with speed vy in the x direction can then be written 



as 



dV V d 

F y = ^^E^HW^ + COmP ' C ° nj -' (9) 

where V projects off the ground state of the vortex. Since dV/dxo is the commutator of H with the partial derivative 
d/dxq, the denominator cancels with the H in the denominator, and so the expression is equal to the Berry phase 
form 

F »— ^asW+^ssW- (10) 

Since the Hamiltonian consists of kinetic energy, a translation invariant interaction between the particles of the 
system, and the interaction with the pinning center, which depends on the difference between the pinning center 
coordinates and the particle coordinates, the derivatives d/dxo, d/dyo, can be replaced by the total particle momentum 
operators — d/dxj, — Y, d/dyj. This gives the force as a commutator of components, P x , P y of the total momentum, 

Fy = -iVv(* To \[P x ,Pyp ro ) . (11) 

At first sight one might think that the two different components of momentum commute, but this depends on boundary 
conditions, since the momentum operators are differential operators. Actually this expression is the integral of a curl, 
and can be evaluated by Stokes' theorem to get 

F 3/ = ^'(* rQ |-i^grad j |* ro >- dr = j> j s ■ dr , (12) 

where the integral is taken over a loop at a large distance from the vortex core. This gives the force in terms of the 
circulation of momentum density (mass current density) at large distances from the vortex. 

Our result that the transverse force is equal to vy times the line integral of the mass current is independent of the 
nature or size of the pinning potential. The general form of this is 

F t = p s K s X vy + p n K n X vy , (13) 

where K„ represents the normal fluid circulation. 

In equilibrium the circulation of the normal fluid around a stationary vortex is zero, since circulation of the normal 
fluid gives rise to viscous dissipation of energy, which in turn leads to growth of the area of the normal fluid vortex 
core. If there is any nonequilibrium normal fluid circulation, it is not obvious that it should be quantized, or that the 
motion of the normal fluid vortex should be correlated with the motion of the superfluid vortex. If only the superfluid 
participates in the circulation round the vortex core, which seems to us to be the most reasonable assumption, this 
gives 

F t = Ps K x vy . (14) 

In combination with the Wexler result for the coefficient of v s , the total transverse force on a vortex is 

F t = p s K s X (vy - v s ) . (15) 

Only the superfluid Magnus force exists unless the normal component participates in the circulation of the superfluid 
to some extent. 

This disagrees with Donnelly's eq. (^) unless the phonon-scattering term proportional to a cancels with the Ior- 
danskii term proportional to p n , and disagrees with Volovik's eq. (|^) unless his coefficient Cp is equal to p n even for 
bosons. The most striking feature of this is that the force is independent of the normal fluid velocity. 

It agrees with the obvious generalization of the classical Magnus force argument to two-fluid dynamics. This 
argument considers the force- momentum balance in a large cylinder surrounding a vortex which is held stationary 
while the fluid flows past it. Bernoulli pressure on the cylinder, and momentum flux across the boundary of the 
cylinder balance with the force on the vortex. In a two-fluid generalization of this there are separate contributions 
from the product of superfluid circulation with superfluid velocity, and from the product of normal fluid velocity with 
normal fluid velocity. 

Since we have only had to consider global properties involving momentum conservation and conditions at a long 
distance from the vortex core, and have not needed to make any detailed consideration of conditions at the core 
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of the vortex, we believe that our arguments are valid for a fermion superfluid described by a single complex order 
parameter. The B phase of 3 He does not quite meet this condition, but the order parameter is an essentially isotropic 
combination of a P-wave orbital state and a triplet spin, so this should behave in much the same way. Volovikll3 ! EJ has 
argued that spectral flow of the unpaired states in the vortex core of-a fermion superfluid leads to a contribution to 
the transverse force that cancels most of the Magnus force, but StoneOJ has examined this argument more closely, and 
does not find that this mechanism is operative unless there is a background to take momentum from these excitations. 
We do not think that there is such a canceling contribution in a homogeneous fermion superfluid. 



V. FORCES DUE TO PHONON SCATTERING 



The result obtained in the previous section, that the coefficient of the vortex velocity in the transverse force is equal 
and opposite to the coefficient of the superfluid velocity leads to the surprising conclusion that normal fluid flow does 
not affect the forca-on the vortex, .unless there is also normal fluid circulation round the vortex. This is surprising, 
because PitaevskiiEJ and IordanskiiEJ argued that the asymmetrical scattering of rotons or phonons by vortices should 
lead to a transverse force when the vortex moves relative to the normal fluid component. 

In the low temperature limit an explicit calculation of the phonon-vortex scattering can be made, and the literature 
quotes a transverse force proportional to p„K x (vy — v„). There are two problems with this result: 

1. The derivation assumes that the phonons interact only with the vortex, but in our argument we assume that 
the phonons, which make up the normal fluid, must be in equilibrium with one another. 

2. In papers from Cleary (1968)0 to Sonin (1997)0 the expression for the transverse force, which is proportional 
to 



T 3 



■(v v - v n ) ^ sin<5 m sin<5 TO+ i sin(£ m+ i - S m ) , (16) 



h 2 <? 

m 



has been rewritten as 



2^3 

(vy - v„) sin 2(5 m+ i - S m ) . (17) 



This would be fine, except that S rn does not tend to zero. If one substitutes the formula 

S m oc sign(m) kqT /Tic 2 , (18) 

which is correct to lowest order in temperature T, into the original formula, a result is obtained which is (at 
least) cubic in Ko and of sixth power in T, or 3/2 power in p n . The second expression is obtained from the first 
by canceling two divergent series, and this gives the quoted expression which is linear in kq and linear in p n \ we 
cannot see any justification for a term of this magnitude. 



VI. SUPERCONDUCTIVITY 



The situation for the transverse force on a vortex in a superconductor is even more confused than the situation for 
a neutral superfluidr-,In the 1960s, Bardeen and StephenO argued for a very small Magnus force, but an analysis by 
Nozieres and VinenBa of an idealized model of a superconductor gave the full value of the Magnus force suggested by 
classical hydrodynamics. 

Wexler's arguments for the coefficient of v s can be applied to the case of a superconductor. When the substrate 
velocity and the vortex velocity zero in the presence of a superfluid electron velocity, this argument gives the expected 
result that there is a Lorentz force on the vortex equal to the integral of ep e v s x B, where p e is the conduction electron 
density. . 

To find the coefficient of vy , Geller, Wexler and ThoulesscB have adapted the arguments of sec. 4 to the very idealized 
model of axharged system with a uniform positive background, which is essentially the situation considered by Nozieres 
and VinencJ, although they, unlike us, also had to assume that the superconductor was extreme type II. This is not 
completely straightforward, even though we have taken the uniform positive background so that we can continue to use 
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momentum conservation, because any choice of the gauge field which is used to describe magnetic effects breaks the 
explicit translation invariance, and makes the implicit translation invariance obscure. Rather than introduce a gauge 
field, we can write the electromagnetic interactions in terms of a Coulomb interaction between electrons and between 
electrons and positive background, together with an instantaneous current-current interaction between the electrons. 
Darwin showed that this is correct up to second order in electron velocity, apart from a relativistic variation of the 
mass with velocity which is unimportant for this problem. We also need a Galilean invariant attractive interaction 
between the charges to produce a paired superconducting state. This gives a Hamiltonian with explicit translation 
invariance, so the arguments of sec. 4 can be taken over. 

The result for the coefficient of vy is formally unchanged, and is the line integral of the canonical momentum 
density on a loop which surrounds the flux line at a distance which is large compared with the penetration length. 
This is actually a surprising result, as at these distances there is no magnetic field or current density produced by the 
vortex line, and the integral is related to the Aharonov-Bohm effect rather than to any classical quantity. 

Since the integral is equal to the trapped magnetic flux, the transverse force can be written as 

p e (v v -v e )xBA, (19) 

so that the transverse force depends only on the motion of the vortex relative to the electrons. It can be rewritten, 
in a form that makes its physical origin more transparent, as 

F f = JJJ p e (v p - v e ) x Bd 3 r + JJJ p e (vv-Vp) xBA. (20) 

The first term is the Lorentz force given by the interaction of the electric current, which is a Galilean invariant, with 
the magnetic field. The second is a Magnus force that acts on the positive substrate moving with velocity v p relative 
to the vortex. The moving vortex generates a dipolar electric charge distribution, which in turn produces a dipolar 
elastic stress on the positive substrate, and this leads to a net force on the positive substrate. 



air results were essentially the same. There have 



A similar analysis was carried out by Nozieres and VinencJ, and tl 
been recent measurements made by Zhu, Brandstrom and SundqvistE3 that support a fairly large value of the Magnus 
force. 



VII. CONCLUSIONS 



We have succeeded in determining the transverse force on a vortex in a neutral superfluid under assumptions that 
are both general and reasonably realistic. Like most other exact results in quantum many-body theory, they are 
related to general conservation laws, and apply, in a slightly different form, to classical systems as well. 

Our generalization to superconductors is far from realistic, since it relies on the uniformity of the positive substrate. 
With a uniform positive substrate an electron gas has infinite conductivity, even in the absence of a pairing interaction, 
so our results can only form a first step towards a plausible theory of the Magnus force in superconductors. We may 
be able to extend the reults from a uniform substrate to an ideal periodic sustrate, but even that is quite inadequate 
for the description of a real metal. We have to be able to take the next step of considering disorder, but there is no 
chance that we will be able to get exact results in that case. 

It would also be interesting to generalize these results to finite systems, nonzero frequency of the vortex motion, 
and a finite density of vortex lines. 

Another line that we are pursuing is the connection between the Magnus force and the quantization of the vortex 
line. We know that there is an intimate connection between the strength of the Lorentz force on an electron and the 
density of degenerate levels of an electron in a magnetic field, and there are good reasons to think that there is a 
similar connection between the strength of the Magnus force on a vortex and the density of degenerate levels for a 
vortex. 
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